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A characteristic lifelength is defined by which a Gaussian distribution is fit to data 
correiated over a 3 sensor array sampling streamwise sidewall pressure. The data were 
acquired at subsonic, transonic and supersonic speeds aboard a Tu-144. Lifelengths are 
estimated using the cross spectrum and are shown to compare favorably with Efimtsov’s 
prediction of correlation space scales. Lifelength distributions are computed in the time/ 
frequency domain using an interval correlation technique on the continuous wavelet 
transform of the original time data. The median values of the lifelength distributions are 
found to be very close to the frequency averaged result. The interval correlation technique is 
shown to allow the retrieval and inspection of the original time data of each event in the 
lifelength distribution, thus providing a means to locate and study the nature of the coherent 
structure in the turbulent boundary layer. The lifelength data can be converted to lifetimes 
using the convection velocity. The lifetime of events in the time/frequency domain are 
displayed in Lifetime Maps. The primary purpose of the paper is to validate these new 
analysis techniques so that they can be used with confidence to further characterize coherent 
structure in the turbulent boundary layer. 


I. Introduction 

The presence of coherence structure in the turbulent boundary layer has been studied for decades, see Robinson 1 
and Gad-el-Hak 2 for excellent overviews of the topic. The promise of coherent structures is that an understanding of 
their behavior would unlock the dynamics of the turbulent boundary layer. While much progress has been made at 
low Reynolds numbers, few results have been obtained at high Reynolds numbers due to difficulties in acquiring data 
in this regime either experimentally or through simulation. The technique introduced here uses wall pressure data and 
an extension of an interval correlation approach introduced earlier 3 to estimate the size and lifetime of the coherent 
structures. The paper deals primarily with validating the approach as some of the results may be considered atypical 
and confidence in the technique must be gained before conclusions are drawn. Most established procedures to which 
this technique can be compared are based on frequency averaged data. The key benefit of the interval correlation 
technique is that the coherent structure can be located and the instantaneous behavior studied in the time domain. 

Coherent structures are located by estimating the characteristic lifelength of the data in the prescribed sub- 
interval of time. The term characteristic lifelength is coined as it is the value of a parameter which fits a Gaussian 
distribution to data correlated over a pressure sensor array. If the analysis where performed in time, rather than space, 
it would be called a lifetime. A related lifetime can be computed from the lifelength using the convection velocity. 
The characteristic lifelength is shown to be related to Efimtsov’s 4 correlation space scales with results comparing 
favorably with Efimtsov’s frequency averaged predictions. The analysis is extended to the time/frequency domain 
using a continuous wavelet transform (CWT). Once in the time/frequency domain, distributions of lifelengths can be 
calculated at any frequency providing additional insight into the coherent structure phenomena. The event lifetimes, 
as derived from the corresponding lifelengths, are plotted in the time/frequency domain and called Lifetime Maps. 
The Lifetime Maps are feature rich and exhibit pronounced activity in concentrated areas and no activity in others, a 
behavior reminiscent of the bursting phenomena described by other researchers 1,2 . 

The discussion begins with a description of the data set. The data were not taken with the subject analysis in mind, 
so the analysis was necessarily adapted to the data. A short array length (3 sensors) is the primary reason for having to 
assume a Gaussian distribution in order to estimate the correlation length. The use of the Gaussian distribution is then 
validated by comparing lifelength distributions derived from frequency averaged cross spectra to Efimtsov’s 
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correlation space scale predictions. The interval correlation technique is introduced with two variants, one using the 
center sensor of the array as a reference, the other using the first principal component as computed from a Singular 
Value Decomposition (SVD) of the data. The reader may be more familiar with the term Proper Orthogonal 
Decomposition (POD) which is an expression common in the literature describing an equivalent computation. 
Typical results obtained in the time/frequency domain are presented and validated by comparing the median 
lifelength frequency distributions to the cross spectrum distributions obtained earlier. The paper concludes with 
example Lifetime Maps and some brief discussion as to what the features observed in the maps may mean. As 
mentioned earlier, the primary focus in this paper is to validate the approach. More work remains to be done to 
adequately explain the behavior displayed in the Lifetime Maps. 

II. Data Set 

The data set was taken aboard a TU-144 as part of NASA’s High Speed Research (HSR) program 5 . Arrays of 
kulites were installed in window blanks positioned in the forward, center and aft sections of the aircraft as shown in 
Fig. 1. The analysis was performed on 3 streamwise kulites in each window blank as indicated in Fig. 1. The sensor 
spacing was 6 cm with the exception that the center-to-upstream sensor spacing was 4 cm on the forward and aft 
window blanks. 
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Figure*l. ‘Sensor layout on Tu-144. Three sensor sets circled 

Three flight conditions were available at Mach numbers 0.6, 0.9 and 2.0. Details of each flight condition are listed 
in Table 1. The freestream velocity, U inf , is computed from the aircraft indicated airspeed. The friction velocity, U T , is 
estimated as 0.03U inf . The convection velocity, U c , was calculated from the delay as indicated by the sensor cross 
correlation. For comparison purposes the common estimate for U c of 0.7U; n f is also given. These values compare very 
well except for the 2.0Ma forward sensor case. At these speeds, one sample delay is equivalent to 80 m/s and thus 
contributes to the large error. The Reynolds number and boundary layer thickness were taken from Aerospaceweb’s 
Atmospheric Properties Calculator which is based on U.S. Standard Atmosphere 1976. Note the high Reynolds 
number in the range of 10 8 for all flight conditions and locations. The boundary layer thickness, 5, ranges from 16 cm 
for the 0.6Ma forward sensor case to 37 cm for the 2.0Ma aft sensor case. 

Results derived from two sets of data will be presented. The first is taken from the forward sensor set at the 3 
Mach numbers. The second, at 0.9Ma, is taken from the 3 different locations. This results in a total of 5 data sets as 
indicated by the rows highlighted in green in Table 1. As discussed above, the convection delay may have 
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considerable error at high Mach numbers. However this is not the only source of error as will be discussed in the next 
section. 


Table 1: Flight Conditions. Selected cases in green. 


Mach 

Altitude 

(km) 

U inf 

(m/s) 

Ut 

(m/s) 

Location 

Distance 

(m) 

U c measured 
(m/s) 

0-7*U inf 

(m/s) 

Re 

5 

(m) 

0.6 

5 

183 

5.5 

Forward 

18.9 

143 

128 

1.6e8 

0.16 

0.6 

5 

183 

5.5 

Center 

32.6 

150 

128 

2.7e8 

0.25 

0.6 

5 

183 

5.5 

Aft 

49.3 

133 

128 

4.1e8 

0.35 

0.9 

9 

273 

8.2 

Forward 

18.9 

200 

191 

1.6e8 

0.16 

0.9 

9 

273 

8.2 

Center 

32.6 

200 

191 

2.8e8 

0.25 

0.9 

9 

273 

8.2 

Aft 

49.3 

200 

191 

4.2e8 

0.35 

2.0 

16.9 

590 

17.7 

Forward 

18.9 

500 

413 

l.le8 

0.17 

2.0 

16.9 

590 

17.7 

Center 

32.6 

400 

413 

2.0e8 

0.27 

2.0 

16.9 

590 

17.7 

Aft 

49.3 

400 

413 

3.0e8 

0.37 


A. Data Quality 

The data were analyzed using an interval correlation technique described below and in Ref. 3. The analysis 
applies an estimated convection delay to the time data in order to emphasize correlated content. This approach 
depends on the assumption that the convection delay is constant both over time and frequency. The former 
assumption is adequate for the purposes of the subject analysis, however, the later assumption has been found to 
cause considerable error under some conditions. 


Correlation between Sensors 2 and 3, Center, 0.9M 
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Figure*2. 'Cross correlation (a) and cross spectrum delay (b) for center section sensors 2 and 3 at 0.9Ma 


Fig. 2 illustrates a best case example. Here the convection delay is computed between the 2 nd and 3 rd center 
section sensors at 0.9Ma using both cross correlation and cross spectrum. The peak delay of 11 to 12 samples in the 
cross correlation (a) is shown to compare favorably with the delay as derived from the phase angle of the cross 
spectrum (b) down to about 300 Hz where the delay begins to oscillate and eventually diverges. This performance is 
adequate as in most cases the behavior of interest is occurring above 300 Hz. 

In contrast, the delay as computed from the same center section at 0.9 Ma but using the 1 st and 2 nd sensors is 
shown to begin to diverge at 700 Hz (Fig. 3b). Even though the cross correlation graph in (a) appears normal, the 


3 


American Institute of Aeronautics and Astronautics 


cross spectrum phase delay in (b) reduces to 10 samples at 1200 Hz and then begins to oscillate at 900 Hz before 
diverging at 700 Hz. Again it should be noted that the data used for sensor 2 in Fig. 3 is identical to the data used for 
sensor 2 in Fig. 2. The cause of this behavior is unknown, but an effect on the analysis will be observed. 


Correlation between Sensors 1 and 2, Center, 0.9M 
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Figure*3. ‘Cross correlation (a) and cross spectrum delay (b) from center sensors 1 and 2 at 0.9Ma 


Phase delay divergence can be found in all the data. For example, the data in Fig. 2(b) may be said to diverge at 
300 Hz. Defining the point of divergence as a cut-off frequency for the analysis, the sensor pairs and related cut-off 
frequencies are listed in Table 2. Four of the fifteen sensor pairs in the 5 selected cases have cut-offs above 500 Hz, 
the worst case being aft sensor pairs 1 and 3 at 2.0 Ma with a cut-off at 1200 Hz. The impact of this behavior will be 
noted where applicable. 

Table 2: Phase cutoff frequencies. Selected cases in green. Noted cases in red 


Mach 

Location 

1>2 

Cutoff 

(Hz) 

2>3 

Cutoff 

(Hz) 

1>3 

Cutoff 

(Hz) 

0.6 

Forward 

300 

300 

300 

0.6 

Center 

800 

800 

700 

0.6 

Aft 

800 

800 

500 

0.9 

Forward 

400 

400 

400 

0.9 

Center 

700 

400 

500 

0.9 

Aft 

400 

400 

400 

2.0 

Forward 

600 

600 

1200 

2.0 

Center 

800 

800 

500 

2.0 

Aft 

300 

400 

600 


4 


American Institute of Aeronautics and Astronautics 


III. Correlation Space Scale 

Efimtsov 4 defines a correlation space scale, A, as the integral of the cross spectrum, <p. in the streamwise direction 
as given in Eq. 1, below. 


A = J |tp(jc, 0, oo)| dx (1) 

o 


Computing this integral over the short distance covered by the three sensors on each blank would not provide a 
meaningful result. Fortunately Efimtsov has derived a general form for A based on the Strouhal number. 


A 


0. 1 Sh \ 2 

u/uj 


+ 


72.8' 


Sh 2 + (72.8/ 1.54) 2 


— 1 / 2 \ 


( 2 ) 


Where 

Sh= to 8U t and, (3) 

U= 0.03 (4) 

Sh is the Strouhal number, U x is the freestream velocity, U c is convection velocity, and, U x is the friction 

velocity. The correlation space scales for the 5 selected cases were computed using the flight parameters from Table 1 
and are shown in Fig. 4. In general the frequency at which the peak in correlation length occurs is seen to increase 
with Mach number, (a), and decrease with streamwise distance, (b). 

Efimtsov Correlation Length Prediction: Forward Location Efimtsov Correlation Length Prediction: 0.9M 




Figure*4.*Correlation space scales for 3 Mach nos. at Forward location, (a), and 3 locations at 0.9 Ma, (b) 


A. Estimated Correlation Space Scale 

As mentioned above, direct computation of the integral in Eq. 3 using the data acquired from the small arrays 
would be impractical. An alternative approach is taken where a characteristic lifelength, X, is defined and optimized 
to fit a Gaussian distribution to the sensor cross spectra referenced to the center sensor. Referencing the center sensor 
ensures a magnitude peak over the center of the array and is key in locating coherent structure when the technique is 
extended to the time domain. The form of the estimated cross spectra distribution is shown in Eq. 5. 
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cp = e 


( 5 ) 


-(x/X ) 2 



Figure'5. 'Center referenced cross spectra (blue) and 
Gaussian fit (red). 


The characteristic lifelength at each frequency is found 
by applying a least squares fit of the Gaussian 
distribution to the cross spectra data. Typical results 
are shown in Fig. 5. 

GO 

A = J cp(jc, X, co )dx (6) 

—00 


A = Xjn (7) 

The integral of the estimated cross spectrum 
distribution, Eq. 6, conveniently evaluates to Xjn , Eq. 

7. The estimated correlation space scale, A , will be 
referred to as the lifelength for the remainder of the 
paper. 


In Fig. 6, the predicted correlation space scale is compared to the measured lifelength for the forward sensor 
cases. The 0.9 Ma case (red) agrees very well with predictions attaining a peak of about 20 cm at 750 Flz vs. the 
predicted 22 cm at 1 kHz. The measured 0.6 Ma (blue) case has a broad peak from 500 to750 Hz and peaking at 14 
cm. This is the correct frequency range, but substantially lower than the 23 cm predicted. The measured 2.0 Ma case 
(green) peaks at 16 cm around 1500 Hz which is substantially lower in frequency and magnitude than the predicted 
case at 25 cm and 2300 Hz. It should be noted that this case was marked in Table 2 as having high cut-off frequencies 
and, therefore, could be expected to return questionable results. Note also that since the 0.9 Ma and 2.0 Ma cases peak 
at lower frequencies than predicted, it might be expected that the 0.6 Ma case would peak at a frequency lower than 
500 Hz very near the data’s cut-off frequency. 


Efimtsov Correlation Length Prediction: Forward Location 
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Figure'6. 'Comparison of predicted correlation space scale, (a), with lifelength, (b), for forward cases 


Fig. 7 compares the 0.9 Ma cases. Here the forward case (blue) is the same as the forward 0.9 Ma case in Fig. 6, 
comparing very well with prediction. The aft case (green) peaks at about the right frequency, but at a much lower 
value, 33 cm vs. 46 cm. Note that the aft case 400 Hz cut-off may have something to do with the lower peak. The 
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center sensor case (red) compares poorly with prediction exhibiting no peak and attaining a maximum value less than 
15 cm versus the nearly 35 cm predicted. Again, the cut-off frequency at 700 Hz for this case may be interfering with 
the analysis as it occurs very near the predicted peak. 

Given the limitations of the data set, the estimated correlation space scale, or lifelength, has been shown to be a 
valid measure. The technique is applied to time domain data using interval correlation as explained in the next 
section. 

Efimtsov Correlation Length Prediction: 0.9M LifeLengths Mach 0.9 




Figure*7. •Comparison of predicted correlation space scale, (a), to lifelength, (b), for Mach number cases. 


IV. Time Domain 

To characterize coherent structure in the turbulent boundary layer, time intervals must be identified in which a 
suspected coherent structure has occurred. Assuming that coherent structures have long lifelengths, this would 
involve scanning the data sets for time intervals in which long lifelengths are detected. Two techniques have been 
developed for this purpose. Both techniques use an interval correlation function, Eq. 8, on data taken from a subset, or 
window, in the time history. One technique uses the center sensor as a reference, the other a value derived from a 
Proper Orthogonal Decomposition (POD) of the data. 


Y} d r*td i 

lC, ' ! 2 , i rms(td t )rms(td 2 ) ^ 

Here the time data, td, is delayed or advanced with respect to the center sensor according to estimates of the 
convection velocity and sensor separation. The summed product of the time data is normalized by the RMS of the 
reference data and the center sensor data. Using the center sensor to normalize the interval correlation was useful in 
rejecting degenerate solutions. This was more of a problem for the POD based analysis as will be described next and 
in section V.E. Comparing POD to Center Referenced. 


A. Proper Orthogonal Decomposition 

Proper Orthogonal Decomposition (POD) is also referred to as Singular Value Decomposition (SVD) and 
Principal Component Analysis. The SVD function is used to break the time delayed data down into its respective 
singular vectors and singular values. The first singular value is then used to construct the first principal component, or 
first mode, of the data. 
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[ t/, *S, V\ = svd(Xtd-y',td 2 ',td^\) 


( 9 ) 


PC, 
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5 ( 1 , 1 ) 
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( 10 ) 


= Z />c i * ,d > 

ic PC u 2,i rms(PC l )rms(td 2 ) ^ ^ 

In Eq. 9, the SVD function is shown in MATLAB notation. The time delayed data, td v are formed into a matrix 
and passed to the SVD function. The values returned are the left and right singular vectors, V and V respectively, and 
the singular values, S. The first principal component, PC b is formed by using only the first singular value S(l,l) with 
the singular vectors as shown in Eq. 10 where the H operator indicates complex transpose. The PC referenced interval 
correlation can then be formed as in Eq. 11. 

An example comparing the first principal component to the data from which it was derived is shown in Fig. 8. 
Note that the SVD does a very good job in constructing a wave shape that has a lot in common with its parent data. A 
very real problem exists with using this kind of decomposition with so few sensors in that correlation is found where 
none really exists. This will be discussed in more detail in section V.E after the time/frequency domain analysis is 
introduced. 


Flow at U,. 


At 2 


At, 



Figure*8.*Time delayed pressure data compared to first principle component. 


B. Gathering Statistics 

Using the interval correlation function on a defined sub-interval of the data set, a number of samples can be 
collected characterizing the lifelength distribution of the data. Fig. 9 shows a typical result for the forward sensor 0.6 
Ma case. The distribution is characterized by its median, maximum and 75 th percentile values. These values for the 
selected cases are shown in Table 3 where the cases with high cut-off frequencies are indicated in red. The results are 
in line with the cross spectrum values, although it is difficult to make a comparison because the time domain result is 
an average over all frequencies. The time domain analysis has another drawback in that the selection of an analysis 
interval skews the results with large intervals masking high frequencies and short intervals biased by very low 
frequencies. This problem is avoided by performing the analysis in the time/frequency domain as will be discussed 
next. 
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frequency, Hz 


Forward, mach: 0.6, corrType: PC.corrLim: 2; window: 100; #Samples: 565 



cm 


Table 3: Time domain lifelength statistics 


Case 

Max 

Median 

75 th pet 

(cm) 

(cm) 

(cm) 

0.6Ma, Forward 

16 

18 

20 

0.9Ma, Forward 

15 

18 

22 

0.9Ma, Center 

17 

18 

21 

0.9Ma, Aft 

19 

21 

25 

2.0Ma, Forward 

14 

18 

21 


Figure*9.*Lifelength distribution taken from 
Forward, 0.6Ma case. 



time, msec 

Figure*10.*Example of CWT of TBL pressure data 


V. Time/Frequency Domain 

The time data can be transformed to the time/ 
frequency domain using a continuous wavelet 
transform. An example transform is shown in Fig. 10. 
More details on this procedure can be found in Ref. 3. 
Once in the time/frequency domain, interval 
correlations can be done at each frequency, i.e., one 
row of data in Fig. 10. Besides the ability to view the 
frequency distribution of the lifelength statistics, this 
approach offers the ability to size the analysis sub- 
interval to a number of periods, thus varying the sub- 
interval size with frequency and removing the bias a 
constant sub-interval introduces. This section is 
largely dedicated to validating the approach, 
identifying its pitfalls and recommending 
improvements. 


A. Lifelength Frequency Distribution 

The lifelength frequency distribution for the forward 
0.6 Ma case is shown in Fig. 11 (a) with the 
distribution’s median and 95 th percentile in (b). In (a) the color in the plots relates to number of events with red being 
highest. It can be seen that as the distribution broadens at lower frequencies, that fewer events occur along the median 
value. A second order effect is that as more events occur at larger lifelengths, fewer events will occur in general. 
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Figure'll. 'Lifelength frequency distribution, (a), and related median and 95 th percentile, (b) for the forward 

sensor 0.6Ma case. 


B. Comparing Interval Correlation to Cross Spectrum Results 

The median lines from the lifelength frequency distributions of the forward sensor cases are plotted in Fig. 12(a) 
alongside the cross spectrum results in (b). Excellent agreement is obtained between these two very different 
techniques. The 0.9 Ma case in (a) peaks at a lower frequency, 400 Hz, than in (b), at 700 Hz, reemphasizing a trend 
away from the Efimtsov prediction showing a peak at 1000 Hz, Fig. 4(a). The 0.6 Ma case shows a similar trend, 
indicating a peak in the median at lower frequencies than predicted. The 2.0 Ma case in (a) follows the cross spectrum 
results very closely. All cases in (a) exhibit an abrupt drop-off below 400 Hz. This is most likely due to errors in the 
phase delay at low frequencies as described earlier. See “Data Quality” on page 3. 




Figure'12. 'Median lifelengths from CWT analysis, (a), compared to lifelength estimates from cross spectrum, 

(b), for the forward sensor cases 

The median of the lifelength frequency distributions for the 0.9 Ma cases are shown in Fig. 13(a) along with the 
related cross spectrum result in (b). One immediately apparent difference in the results is seen for the aft sensor case 
where the interval correlation result does not return that higher lifelength below 400 Hz as indicated by both the cross 
spectrum results in Fig. 13(b) and in Fig. 4(b). The error again seems to be related to the departure of the phase delay 
at the lowest frequencies from the practically constant value observed at mid to high frequencies. 
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1. Resolving the low frequency error 

A simple fix to the analysis would be to introduce a frequency dependent delay to compensate for the observed 
deviation at low frequencies. However, this would be done in the absence of a physical explanation for the observed 
behavior. One possible explanation for the observed phase is that there is a rotation in the flow, causing added phase 
shift as the structure propagates. This would support both a constant convection velocity and a variable phase delay. 

Another possible explanation is that the low frequency behavior is caused by the injection of low pressure, 
freestream fluid into the inner layer. This phenomena is described as part of the bursting process by many who have 
studied coherent structure 1,2 and will not be detailed here. If indeed, freestream fluid is impinging on the inner layer, 
it would be of lower velocity and, most likely, lower frequency, thus providing an explanation for the observed 
behavior. 



LifeLengths Mach 0.9 



200 400 600 800 1000 1200 1400 1600 1800 2000 

frequency, Hz 

Figure*13. 'Median lifelengths from CWT analysis, (a), compared to lifelength estimates from cross 

spectrum, (b), for 0.9 Ma cases. 


C. Comparing Wavelet Transform Data to Time Domain Data 

Once a sub-interval with large lifelength has been located, the original time data can be retrieved. For the 
purposes here, this is done to increase confidence in the analysis. In future analyses this can be done to increase 
understanding of the physical phenomena. 


Start Index: 11521; Frequency: 331.1558 LifeLength: 0.75835 
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Figure*14.'CWT, (a) and wavelet data, (b), taken from selected time interval at 331 Hz. 
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The lifelength frequency distributions are composed of discrete events. Each dot in Fig. 1 1(a) can be queried and 
the originating data displayed. For example, the CWT in Fig. 14(a) is associated with an event detected at 331 Flz 
with a lifelength of 0.75m. The CWT is 3 sub-intervals long with the subject data highlighted in the figure. Several 
features are apparent in the data. One is that there appears to be a continuous tone just below 400 Hz which morphs 
into the subject event. Another feature is the 2 nd harmonic at ~660Hz which appears in the center sub-interval. One 
might wish to know the lifelengths of the features surrounding the detected event. This is the subject of the following 
section, Lifetime Maps. The actual wavelet filtered data is shown in Fig. 14(b). The high degree correlation in the 3 
sensors which led to the long lifetime is obvious. To test whether this is real or somehow a result of the analysis we 
look to the original time data. 

The time data associated with the graphs in Fig. 14 is shown in Fig. 15(a). The presence of a persistent low 
frequency can be seen in the data with a more pronounced feature centered in the detected sub-interval. 


Time History Time History. 200 Hz BP. Power Fraction: 0.33 



Figure*15. ‘Original time history data, (a), and result of 200 Hz bandpass filter, (b). 

The bandpass filtered time data is shown in Fig. 15(b). The filter is centered at 331 Hz with a bandwidth of 200 
Hz. This is wide compared to the 10 Hz binwidth of the CWT, but similarities can be seen in the data, i.e., comparing 
Fig. 15(b) with Fig. 14(b). Notice, for example, that between 4 and 8 msec, the red curve is larger and leading, and, 
the green curve is smaller and lagging in both figures. Also note the between 16 and 18 msec that the amplitude and 
phase of the 3 signals are almost identical in the two figures. These kind of comparisons improve with narrower filter 
bandwidth and have been found to be consistent throughout the test cases. This raises confidence in the technique to a 
high level encouraging further analysis and insight. For example, note that the bandpass filtered data is shown to have 
33% of the RMS power contained in the original data. One percent of the bandwidth contains 33% of the power. This 
indicates that coherent structures do contain a lot of power relative to less correlated phenomena. 

D. Event Detection Rate 

An event detection rate can be obtained by the ratio of the total number of events detected to the total number of 
intervals sampled. This figure is not an estimate of how often an event occurs, but how often it is detected. The 
difference here is that if the sample interval is significantly smaller than the feature size, the event will be detected 
several times, thus increasing the detection rate. This is discussed in more detail the following section. Lifetime 
Maps. The detection rates for the 0.9 Ma cases are shown in Fig. 16(a) with the corresponding lifelength median and 
95 th percentile in Fig. 16(b). The center section has the highest rate and smallest lifelength while the aft section has 
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the lowest rate and largest lifelength. This relationship is as expected in that longer lasting events will occur less 
often. 




Figure*16. 'Detection rate, (a), and, lifelength median and 95th percentiles for 0.9 Ma cases. 


E. Comparing POD to Center Referenced. 

The data presented above are largely from interval analysis using PC! referenced interval correlation. The POD 
decomposition has the ability to generate a wave shape with the highest correlation among the data sets. The Center 
Referenced analysis works with the data as they are and thus consistently returns lower lifelength values. This effect 
can be seen in Fig. 17(a) where the POD and Center Referenced (CR) analyses are compared for the forward sensor 



Figure»17.'POD and Center Referenced results compared. Forward sensor 0.6 Ma (a), random data (b). 

0.6 Ma case. Note that the two results are very similar with the exception that the CR results are lower. This brings up 
the question as to whether the indicated coherent structure is real or serendipitous. Are actual long lifelength features 
being detected or are they just the accidental overlay of otherwise random events. To test this, both analyses were run 
on completely random data. The results, Fig. 17(b), indicate that the POD will indeed extract coherence from nothing. 
The Center Referenced approach is less susceptible on average, but in the extreme, i.e., for long lifelengths, both will 
detect events where there should be none. One reason for this is that this particular data set consists of only 3 sensors. 
Increasing the array size both in number and dimension will increase the ability of the analysis to reject random 
events. 
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VI. Lifetime Maps 

The lifelength data derived from the interval correlation analysis can be converted to lifetimes and plotted in the 
time/frequency domain. The lifelengths are converted using the estimated convection velocity. The resulting Lifetime 
Maps can then be conveniently compared to the continuous waveform transform plots. Examples of these plots are 
shown in Fig. 18 through Fig. 20 for the 0.9 Ma cases. In the Lifetime Maps, the color of the indicator is the lifelength 
and the length of the indicator in the time direction is the lifetime. The analysis intervals are 4 cycles long and are 
sampled every 2 msec so that the samples are continuous at 2000 Liz and overlapped 50% at 1000 Liz. The Lifetime 
Map, Fig. 18(b) appears to have lines of constant time at high frequencies which bend towards increasing time at low 
frequencies. This is due to the lifetime indicators being plotted in the center of the analysis interval which increases in 
size at lower frequencies. This plotting technique was used so that the lifetime indicators would be placed close to the 
CWT data from which it was derived. Lifetime data derived from the center sensor are presented first as it displays 
clear indicators of behavior which become less obvious in the forward and aft cases. 



Figure*18.*CWT, (a), and lifetime map, (b), taken from center sensors at 0.9Ma. Color in (b) relates to 

lifelength in meters. 

The dominant feature of the Lifetime Map of the center sensor in Fig. 18(b) is the area from 25 to 50 msec and 
400 to 600 Liz. This area contains a region of lifetimes which approach the sample period of 2 msec, implying 
continuous correlated activity. The long lifetime region corresponds to a feature in the CWT at 35 msec and from 400 
to 600 Liz which resembles a flame, having a broad base at low frequency and necking down at higher frequencies, 
becoming almost wispy above 1 000 Hz. These types of features are common in the CWT, but do not always produce 
long lifetime indications. The regions to the left and right of the long lifetime region in Fig. 18(b) are best described 
as two ellipsoid regions void of correlated activity but surrounded by long lifetime events which delineate the voided 
regions. It is common to observe the widening and abrupt termination of a lifetime line as it is traced in frequency. 
Quite often the line reappears after small gap, resuming a long lifetime indication before narrowing down again. 
Several examples of this kind of behavior can be found in all Lifetime Maps. The behavior appears to indicate that the 
data in the voided region are nonlinear, varying in frequency and/or phase, defeating the interval correlation analysis. 

The CWT and Lifetime Map for a sample of data taken from the forward section are shown in Fig. 19. An 
immediate observation is that the Lifetime Map of the forward section is less dense than that of the center section. 
This follows from Fig. 16 which indicates that the center section had the highest detection rate and lowest lifelengths. 
The majority of the behavior in the center section data appears average, i.e., most of the lifetimes are < 1 msec with 
occasional bursts of longer lifetime behavior. The Lifetime Map of the forward section data, Fig. 19(b), contains 
larger voided areas than the center section. This would contribute to a lower detection rate, again in agreement with 
the result plotted in Fig. 16. As in the center section case, the voided areas are surrounded by long lifetime events 
giving the impression that the analysis is failing to detect long lifetime structures that appear to exist in these regions. 
A common characteristic found in voided areas can be observed in the CWT in Fig. 19(a) in the area which 
corresponds to the voided region from 40 to 75 msecs around 400 Hz in the Lifetime Map in Fig. 19(b). What appears 
to be strong 1 st and 2 nd harmonics exist at the borders of the region, where long lifetimes exist, with the interior, 
voided, region becoming bifurcated. It is possible that the void is a normal representation of strong bimodal behavior. 
This kind of modal behavior is in contrast to the long, continuous lifetime lines in Fig. 18(b) which are more 
indicative of broadband behavior with short lifelengths. It is then the voided areas, including the surrounding long 
lifetime borders which might indicate the presence of a coherent structure. The voided regions are at least 10 msec in 
duration. This would correspond to a feature length of 2 m, very much larger than the estimated lifelength of about 
0.2 m. 
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Figure*19.*CWT, (a), and lifetime map, (b), taken from forward sensors at 0.9Ma. Color in (b) relates to 

lifelength in meters. 
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The CWT and Lifetime Map taken from the aft section are shown in Fig. 20. The Lifetime Map follows the 
behavior inferred from Fig. 16, i.e., there are fewer events detected and a greater percentage are long lifetime than 
either the center or forward sections. The aft section appears to be less able to support broadband behavior with most 
lifetime lines interrupted by voids. Below 600 Hz the voided areas appear to be about the same size as those in the 
forward section. At 800 Hz and above, the Lifetime Map of the aft section data becomes sparse, with large voided 
areas. The differences in high frequency behavior as indicated in the Lifetime Maps of these 3 cases may have 
resulted from the apparent differences in energy at the higher frequencies as seen in the respective CWTs. The center 
section appears to have more high amplitude activity above 600 Hz than either the forward or aft section. 
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Figure*20.*CWT, (a), and lifetime map, (b), taken from Aft sensors at 0.9Ma. Color in (b) relates to 

lifelength in meters. 


VII. Conclusion 

Interval correlation analysis has been introduced using both a Center Referenced and first Principal Component 
referenced approach. The analyses calculate the correlation length, or lifelength, of coherent structure in the turbulent 
boundary layer based on an assumption of Gaussian distribution over the sensors. This assumption was shown to be 
valid by comparison to predictions based on Efimtsov’s correlation space scale. The analysis was extended and 
verified in the time/frequency domain using continuous wavelet transforms. The computation produces lifelength 
distributions at each analysis frequency and allows the retrieval of the time data that produced a specific result, thus 
yielding insight into the physical process that averaged statistics cannot provide. It is found that the use of the first 
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Principal Component as reference often results in over estimation of the lifelength, especially in the case where too 
few sensors are available. 

The analysis was extended to include Lifetime Maps where event lifetimes, as computed from the lifelengths, are 
plotted in the time/frequency domain. These plots show evidence of large nonlinear structures, ~ 2 m in length, as 
indicated by regions delineated by long lifetime events. 
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